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ABSTRACT 

Context. The existence of debris disks around old main sequence stars is usually explained by continuous replenishment of small dust 
grains through collisions from a reservoir of larger objects. 

Aims. We present photometric data of debris disks around HEP 103389 (HD 199260), HIP 107350 (HNPeg, HD 206860), and 
HIP 1 1 4948 (HD 219482), obtained in the context of our Herschel Open Time Key Program DUNES (DUst around NEarby Stars). 
Methods. We used Herschel/PACS to detect the thermal emission of the three debris disks with a 3<r sensitivity of a mJy at 
100 /rm and 160/<m. In addition, we obtained HerschellPACS photometric data at 70 pm for HIP 103389. These observations are 
conqrleinented by a large vari^y of optical to far-infrared photometric data. TVo difiTarent approaches are applied to reduce the 
Herschel data to investigate the impact of data reduction on photometry. We fit analytical models to the available spectral energy 
distribution (SED) data using the fitting method of simulated thermal annealing as well as a classical grid search methr^. 

Results. The SEDs of the three disks potentially exhibit an unusually steep decrease at wavelengths > 70/irm. We investigate the 
significance of the peculiar shape of these SEDs and the impact on mc^els of the disks provided it is real. Using grain compositions 
that have been applied successfully for modeling of many other debris disks, our modeling reveals that such a steep decrease of the 
SEDs in the long wavelength regime is inconsistent with a power-law exponent of the grain size distribution -3.5 expected from a 
standard equilibrium collisional cascade. In contrast, a steep grain size distribution or, alternatively an upper grain size in the range 
of few tens of micrometers are implied. This suggests that a very distinct range of grain sizes would dominate the thermal emission 
of such disks. However, we demonstrate that the understanding of the data of faint sources obtained with Herschel is still incomplete 
and that the significance of our results depends on the version of the data reduction pipeline used. 

Conclusions. A new mechanism to [nrxluce the dust in the presented debris disks, deviations from the conditions required for a stan- 
dard equilibrium collisional cascade (grain size exponent of -3.5), antVor significantly different dust properties would be necessary 
to explain the potentially steep SED shtpe of the three debris disks presented. 

Key words. Stars: circumsteUar matter - Stars: individual: HIP 103389, HIP 107350, HIP 114948 - Infrared; planetary systems - 
Infrared: stars 


1 . Introduction 

* Herschel is an ESA space observateny with science instruments 

provided by European-led Principal Investigator conscotia and with im- Debris disks were first discovered via infrared excess emissiem 
potta.-it participation from NASA. associated with main sequence stars such as Vega detected by 
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the Infrared Astronomical Satellite (IRAS; Aumann et al. 1984). 
The first spatially resolved image of a debris disk was that of 
ySPictoris (Smith & Terrile 1984) in optical scattered light. In 
the last few years, the Spitzer Space Telescope has revealed 
that debris disks are common around main sequence stars (e.g.. 
Trilling et al. 2008). Several Key Programs on the Herschel 
Space Observatory (Pilbratt et al. 2010) are dedicated to the 
study of various aspects of the formation and evolution of plane- 
tary systems and their attendant circumstellar debris disks (e.g., 
Augereau et al. 2008). Our Herschel Open Time Key Program 
(OIKP) DUNES (DUst around NBarby Stars; Biroa et al. 2010) 
aims to detect debris disks with fractional luminosities simi- 
lar to the Edgeworth-Kuiper Belt level {LijLi, - 10~^ to 10"®; 
Stem 1996; Vitense et al. 2010) around a volume limited sample 
id < 20 pc) of Sun-like stars (F, G, and K spectral type). Some 
additional sources at20pc<d<25pcare included because of 
their known excesses that have previously been detected with 
Spitzer, or because they are known exoplanet host stars. 

As the most readily detectable signposts of other planetary 
systems, debris disks help us to improve our understanding of 
the formation and evolution of them as well as of our own solar 
system (Meyer et al. 2007; W^att 2008; Krivov 2010). Studying 
the spectral energy distribution (SED) of the dust alone usually 
protides only weak, ambiguous constraints to their properties 
such as chemical composition, grain size, and spatial dikribu- 
tion (Wolf & Hillenbrand 2(X)3). For example, the locaticHi of 
the inner disk radius of the dust distribution is strongly degener- 
ate with the lower limit of the grain size distribution. However, 
the grain size distribution can be described by a power-law and 
the power-law exponent derived from pure SBD fitting is usu- 
ally well constrained and consistent with the analytical value of 
-3.5 derived by Dohnanyi (1969) under several assumptions: (1) 
The particles are produced through a collisional cascade from 
infinitely large to extremely small grains. (2) The grains arc not 
affected by any other influence such as stellar radiation and drag 
fcacss. (3) The strength of the particles, i.e., the energy per vol- 
ume fraction necessary to disrupt them, is independent fom the 
grain size. This is further referred to as a standard equilibrium 
cotbsional cascade. 

In this paper, the potential HerschelfDUNBS discov- 
ery of an unusually steep decrease of the SEDs of three 
spatially unresolved debris disks around Sun-like stars, 
HIP 103389 (HD 199260), HIP 107350 (HNPeg, HD 206860), 
and HIP 114^ (HD 219482), is jaesent^'. It is demonstrate 
on the example of the three targets that the analysis of data of 
faint point sources obtained with Herschel depends very much 
on the photometric calibration and the exact determination of 
the uncertainties. In particular, the significance of our results 
depends on the version of the data reduction pipeline used. 
Provided it is real, the steep decrease occurs in Ae range of 
10 pin to 160/rm, inconsistent with a dust grain size distribution 
follcwing a power-law derived from a standard equilibrium col- 
lisional casc^e. The 100 /cm and 160/rm data represent the first 
photometric measurements of these faint disks at wavelengths 
> 70 ^m and, thus, are the first ones to reveal the unusual shape 
of the SEDs. Results from detailed SED modeling including ad- 
ditional photometric data from the literature are presented as 
well. 


* Two of the sources, HIP 103389 and HIP 107350, are shared targets 
betwsen the DUNES survey in the context of which the analysis is car- 
ried out and the DEBRIS survey (Matthews et al. 2010; l%illips et al. 
20101. 





Fig. 1. Images of the three sources at 100/rm (top) and 160 //m 
(bottom) from Reduction 1 (Sect. 2). The images arc displayed 
in a linear stretch from -0.1 fpeak to of the source of 
interest The peak flux is 0.38 mJy/pixel, 0.23 mJy/pixel, and 
0.54 miy/pixel at 100//m and 0. 17 inJy /pixel, 0.16inJy/pixeI, 
and 0.45 i^y /pixel at 160 /rra for HIP 103389, HIP 107350, and 
HIP 1 14948, respectively. The white circles in the 160^ im- 
ages are centered on the brightest pixel of the source, respec- 
tively. 


Observations and data reduction are described in Sect. 2 to- 
gether with the basic observational results. A theoretical discus- 
sion on how steep the decrease of the SED of a debris disk to- 
ward longer wavelengths is expected to be and a characteriza- 
tion of unusually steep SED sources are given in Sect 3. A de- 
tailed description of our modeling of the systems can be found 
in Sect 4 and results are discussed in Sect. 5. Conclusions are 
given in Sect. 6. 


2. Observations & data reduction 

Tvo Herschel/PACS mini-scan map observations of each target 
were taken with the 100/160 channel combination at array ori- 
entation angles of 70° and 110° providing scan and cross-scan 
coverage to assist in the removal of noise artifacts from the fi- 
nal composite mosaic, fri addition, two scan map observations 
of HIP 103389 were taken with the 70/160 channel combina- 
tion with the same array orientations. Each scan map consists 
of 10 legs of 3' length, with a 4" separation between legs, at the 
medium slew speed (20" per second). In this way, a region of 
s» 1 square arc minute around the source position was covered 
to uniform depth in the resulting mosaic. A summary of the ob- 
serving log is presented in Thble 1. 

PACS data reduction was carried out in version 4.2 
(Reduction I) and 7.2 (Reduction 2) of HIPE (Ott et al. 2010) 
starting from the level 0 products using a modified version of 
the standard reduction script (provided within HIPE). Two dif- 
ferent reductions have been carried out, because changes in the 
still developing pipeline of HIPE may significantly affect the re- 
sults. The separate scans at the two position angles of each chan- 
nel pair were mosaicked to produce a final image at each wave- 
length. Output scales for the final mosaics are 1" per pixel for the 
70/imand lOO/rm images and 2" per pixel for the l^pm image 
(smaller than the native detector pixel sizes of 3'.'2 at 10pm and 
100 pm and 6’.'4 at 160pm, respectively). 

For the further analysis of the data, tfifferent approaches have 
been used for Reductions 1 and 2. These are described in the 
following: 
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Table 1. Summary of ffewcAeZ/PACS observations of the three 
objects 


Object 

Observation ID 

OD 

4[/im] 

T[s] 

HU 103389 

1342193157/58 

322 

70/160 

180 

HIP 103389 

1342193159/60 

322 

lOQ/160 

1440 

mP 107350 

1342195779/80 

355 

100/160 

1440 

HIF 114948 

1342196803/04 

372 

100/160 

1440 


Notes. OD is the observing date that counts the elapsed time in days 
since launch on 14 May 2009, A gives the central wavelength of the 
two filters used simultaneously, and T is the total on-source integration 
time. 

Reduction 1; A high-pass filter was used to remove large scale 
background variation from the images, with filter widths of 31 
frames at 70/rm and 100 ;rm and SI frames at 160^ (equiva- 
lent to ~ 60" and 100", respectively). A central region of 30" 
radius in the images was masked from the high pass filter pro- 
cess to prevent the removal of any faint, extended structure near 
the source position. PACS fluxes and sky noise were measured 
by aperture photometry using a custom script based on the IDL 
APBR routine, batted on the DAOPHOT aperture photometry 
module and cross-checked with the internal HIPE aperture pho- 
tometry routines. The aperture r^us and sky annulus dimen- 
sions were 20" and 30" - 40", respectively at all bands. Sky 
noise for each wavelength was calculated from the rms vari- 
ance of the sky annulus pixel values multiplied by the beam 
size. Results were scaled by the appropriate aperture correction. 
Absolute uncertainties are 5% at 70/im and 100 jum and 10% at 
160/im using the calibration files provided along with this HIPE 
version. 

The images of the three sources from Reduction 1 at 100 /nn 
and 160 /on are shown in Fig. 1. 

Reduction 2: The same filter width as in Reduction 1 has been 
used for the high-pass filtering, but sources were masked based 
on a threshold criterion, the value of which was calculated from 
the standard deviation of all non-zero pixels in the image. PACS 
fluxes and sky noise were measured by aperture photometry us- 
ing both the intcnml HIPE aperture photometry routines and 
the MIDAS data analysis package to check for consistency. The 
source flux was measured within a circular aperture of 4", 5", 
and 8" radius at 70 /rm, 100 /rm, and 160 //m, respectively, and 
scaled by the appropriate aperture correction. Sky noise for each 
wavelength was calculated from the rms variance of ten sky aper- 
tures of the same size as the source aperture and randomly dis- 
tributed across the uniformly covered part of the image. The re- 
sulting deviation of these apertures was calibrated by the aper- 
ture correction and correlate noise correction factors. Extensive 
testing has shown that this approach efficiently accounts for the 
correlated noise in the output images. In fact, this approach pro- 
vides slightly conservative uncertainties. An extended discus- 
sion will be included in the survey overview paper (Biroa et al., 
in prq).). Aperture correction factors have also been qrplied to 
the rms aperture flux dispersion. Photometric calibration uncer- 
tainties in this version of HIPE are 3% - 5%^. 


of HIP 103389, HIP 107350, and HIP 114948 (Eiroa et al., in 
prep.). The stellar contribution to the total SED of each star 
is estimated using a PHOENIX/GAIA synthetic stellar model 
(Brott & Hauschildt 2005). Specific models for Tes, log g* and 
[Fe/H] of each star are built by interpolation from the available 
grid. A normalization is applied to the short wavelength range 
(A < 20 /rm) of the long wavelength section of the Spitzer/IRS 
spectrum where available. Otherwise, the normalization is done 
to the flux measurements at photometric bands from B to K. 
All three sources arc found to have significant excess 3 c) 
at 70/im and 100 /rm. HIP 1 14948 also has significant excess at 
160/rm. HIP 103389 and HIP 107350 have estimated excesses 
at 160/rm between 0.8tr and 2.4<r, depending on the reduc- 
tion. Offsets between optical positions and positions in the PACS 
70/rm images arc below 1 cr of the Herschel pointing accuracy. 
In each case, the observed excess is attributed to the presence ^ 
a debris disk associated with the star. The measured FWHM of 
the sources (derived from a 2-D Gaussian fit) are consistent with 
unresolved objects. One can therefore constrain the extension of 
the emitting area (adopting the distance of the objects listed in 
Table 3) to a diameter of less than 7.' 1 (156 AU) for HIP 103389, 
7.' 1 (126 AU) for HIP107350. and 7'2 (148 AU) for HIP 1 14948. 

The spectral slopes of the three systems form Reduction 1 
and Reduction 2 are listed in Thble 2. The SEDs of all three 
sources are found to exhibit an unusually steep decrease at the 
wavelength range of 70/rm to 160/rm (steeper than a black 
body radiator in the Ray leigh-J cans regime) from the results of 
Reduction 1 . The steepness is evaluated based on the spectral in- 
dex of the SEDs between two different wavelengths (see Sect. 3 
for a detailed description). The significance of this steepness 
is evaluated using error propagation. However, the results from 
the new (but not necessarily better) version of the data reduc- 
tion (Reduction 2) give a shallower decrease of the SED toward 
longer wavelengths. Thus, we emphasize that this is a potential 
discovery of a steep decrease of the SEDs of the three sources 
and that a deeper understanding of the data (in particular for very 
faint sources) than possible until now is necessary to make a fi- 
nal statement. In the subsequent sections, we discuss the conse- 
quences of this steepness of the SEDs under the assumption that 
it is real. 


3. What is an unusually steep SED? 

Infrared excesses observed from debris disk systems are the re- 
sult of absorption of incident starlight by dust grains in the disk 
and subsequent re-emission of stellar r^ation. Since the dust 
is colder than the star, its emission peaks at longer wavelengths 
(we assume that there is only one peak). In this section, a quali- 
tative discussion is given of how steep the SED of a debris disk is 
expected to be. We first investigate the shape of the excess treat- 
ing both the star and the disk as single temperature black body 
radiators. Later, we will discuss the efiects of a more realistic 
disk model, i.e., allowing for a range of temperatures and more 
realistic grain properties. 


2.1. ObservaBonal results 


3.1. Treating the dust as a single temperature black body 


The data obtained in this work for the three sources using 
both reductions as well as the predicted photospheric fluxes 
are listed in Thble 2. Table 3 lists the main stelW parameters 

^ Release note PICC-ME-TN-037 


The spectral index A of an SED is defined as 

^ _ dlOgFy 
dlogv 


( 1 ) 
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servational results obtained in this woric” 

Source 


HIP 107350 

HIP 114948 



Reduction 1 Reduction 2 

Redaction 1 Reduction 2 

FrMxn [mJy] 

47.4 ± 2.7 44.0 ± 2.3 

•ft ... 

... ... 

Fmipsvo [mJy] 

46.6 ±3.8 

28.4 ±2.5 

68.7 ± 3.0 

J^PACSioo [mJy] 

23.7 ±1.4 26.3 ±1.7 

11.0 ±0.9 15.1 ±1.3 

42.5 ± 2.2 40.8 ±1.6 

[mJy] 

5.0 ±1.3 7.7 ±2.5 

4.4 ±1.5 4.4 ±2.3 

12.7 ±1.9 13.3 ± 2.2 

f *, 70 ,/m'’ [mJy] 

13.8 

13.4 

15.0 

Ft,,mm [mJy] 

6.8 

6.6 

7.3 

[mJy] 

2.7 

2.6 

2.9 

Aioo.70 

1.94 ± 0.32 1.44 ±0.33 

2.66 ± 0.45 1.77 ± 0.49 

1.35 ± 0.26 1.46 ±0.23 

Aiso.ioo 

3.31 ±0.69 2.61 ±0.82 

1.95 ±0.88 2.62 ±1.30 

2.57 ± 0.44 2.38 ± 0.43 

Aim,™ 

2.72 ±0.39 2.11 ±0.46 

2.25 ± 0.53 2.26 ±0.74 

2.04 ±0.22 1.99 ±0.25 



7.'1 Xft'O 

7.'2 X t'i 

RAp»csioo 

20“56"47'.47 

21“44“31'.27 

23*’16”57.46 

DECpAcsioo 

-26"17'45'.'5 

+14M6'19'/4 

-62°00'05'.'7 

RAgoUcai 

20'‘56”4r.33 

21'‘44”31*.33 

23'-16"57.69 


-26’’17'47.'0 

+14°46'1970 

-62°00'04<.'3 


Notes. Uncertainties of the flux measurements are total uncertainties including sky noise and calibration uncertainties. A,,,,, is the spectral index 
measured from vi and vj (Sect 3), identified with the corresponding PACS bands. 

The 5pirzer/MIPS 70pm photometry is also listed, since these data complement the Herschel observations and ate important to illustrate 
the unusual behavior of the SEDs. These data ate published by Beichman et al. (2(X)6, HIP 103389 and HIP 1 14948) and Bryden et al. (2006, 
HIP 107350) and have been re-reduced in the context of DUNES (Eiroa et al., in prep.). ® Predicted stellar photosphere using a PHOENIX/GAIA 
synthetic stellar model (Eiroa et al., in prep.). FWHM as measured iioro the PACS images at 100pm using a 2-D Gaussian iiL 


where Fy is the total flux of the object at the frequency v. It is 
typically measured as the slope of the SED between two data 
points available, i.e., 


logfva-log^vi 

log V2- log VI 


( 2 ) 


We further denote the two frequencies vi and V 2 with their cor- 
responding PACS wave bands (e.g., Aiqo, 7 o). The spectral index 
of a black body radiator amounts to A = 2 in the Rayleigh-Jeans 
regime. It is decreasing toward shorter wavelengths where the 
Rayiei^-Jeans approximation is not valid. A spectral index of 
A > 2 means that the SED falls off toward longer wavelengths 
steeper than a single temperature black body in its Rayleigh- 
Jeans regime. In the wavelength range where dust rcemission 
is observed, the emission of the star can be very well approxi- 
mated by a black body radiator in its Rayleigh-Jeans regime. It 
will only contribute as an additional component with a spectral 
index of A = 2. Thus, the contributiem of the star will - if sig- 
nificant - move the spectral slope of the whole system closer to 
2 . 


3.2. More exact treatment of the disk 

Our first step to expand the above discussion is to allow for dif- 
ferent dust temperatures, while the dust is still treated as a num- 
ber of single temperature black bodies (i.e., a radially extended 
disk). A two component black body has been used successfully 
to model SEDs of debris disks (e.g., Hillenbrand et al. 2008). 
This results in a flatter SED (Avi,,^ < 2) than in the above case, 
since the colder dust component will contribute additional flux 
in the long wavelength regime of the warmer dust component. 
Thus, using several components with different temperatures one 
can only reach Ay,,v 2 = 2 in the long wavelength regime of the 
coldest dust component (i.e., at even longer wavelengths com- 
pared to the single black body case above). At shorter wave- 
lengths A„,^ will always be < 2. 



X [#tm] 


Fig. 2. Absorption efficiency QtbtfX) for astronomical silicate 
grains of different grain radii a (see Sect. 3.2 for details). 


Using Mie theory, one can derive absorption efficiencies 
Qoit ('i) for more realistic dust grains. The emission of a dust 
grain is then described by the following equation which is a mod- 
ification of Planck’s law: 


2/ic^ 


exp[hc//Ur*«)-l 


( 3 ) 


The shape of the simple black body SED is modified by the 
function <2ats (<!) which will result in a steeper decrease with 
wavelength in the Rayleigh-Jeans regime of the black body if 
Gtb» (-1) is decreasing with wavelength. Fig. 2 shows gabs (-1) fw 
astronomical silicate (Draine 2003) and different grain radii a. 
f2at» (<1) is close to 1 at short wavelengths, while it indeed ex- 
hibits a break at A » 2jra and then decreases toward longer wave- 
lengths with (or v^). This is analog to the discussion of the 
opacity index, e.g., by Draine (2006). An SED of a dust disk that 
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Tkble 3. Physical properties of the three stars considered in this work (for details and references sec Sect. 4.1) 


Object 

HIP 103389 

HIP 107350 


Name 

HR 8013, HD 199260 

HN Peg, HD 206860 

HR 8843, HD 219482 

Distance (pc) 

21.97 

17.88 

20 J4 

Spectral type and lumioosity class 

F7V 

GOV 

F7 V 

Range of spectral types 

F6V-F8V 

GOIV-V 

F6V-F8V 

V,B-V 

5.70, 0.51 

5.96, 0.59 

5.64, 0.52 

Absolute magnitude Mv, bolometric correction 

3.99, -0.02 

4.70, -0.05 

4.08, -0.02 

Bolometric luminosity, L, [Lo] 

2.025 

1.090 

1.867 

Effective temperature (K) 

6257 

5952 

6240 

Surface gravity, logg* 

4.36 

4.44 

4.31 

Radius, R* [Rol 

1.2 

0.99 

1.17 

Metallicity [PeAl] 

-0.14 

-0.07 

-0.21 

Radial velocity (kn^s) 

-16.1 

18.1 

0.47 

Rotaticmal velocity, vsini (km/s) 

13.7 

12.8 

9.0 

Rotation period, P (days) 

4.4 (Pf sin 0 

4.7 

- 

Li I equivalent width (mA) 

80.0 

91.9 

73.7 

Space velocities U, V, W (kni/s) 

-17.2,-10.2,1.1 

-13.9,-20.1,-11.1 

14.9, -8.0, -5.7 

Activity, logR|g( 

-4.402 

-4.48 

^.434 

X-Ray luminosity, log Lx/Lt 

-4.71 

-4.37 

-4.43 

Mass, [AIq] 

1.28, 

0.98 

1.02 

Age (Gyr) 

0.304, 0.184, 0.748(sin0 

0.530,0.162,0.291 

0.385,0.115,- 


Notes. The ages given are derived using the following methods (left to right): log/{^ activity index levels, ROSAT X-Ray himinosities /.x/f^oi. 
Rotation pmod. For more details and references, see Sect 4.1. 


is steeper in the long wavelength regime than a black body is ob- 
viously possible. However, debris disks are usually expected to 
be radially extended (e.g., Backman et al. 2009; Ertel et al. 201 1 ; 
Lohne et al. 2011) and to contain particles of different size. This 
suggests a broad range of dust temperatures to be present result- 
ing in a flattening of the SED. In addition, the grain size depen- 
dence of the break in {2«i» W means that the grains have to be 
smaller than » Allut to result in an SED that falls off steeper than 
a black body (Draine 2006). This suggests grains smaller than 
~ 15/nn for the three disks presented in this paper (Table 2). 
Thus, one qualitatively expects a narrow dust belt composed of 
small grains to be present. While a spectral slope steeper than 
2 is common for ^bris disks at (sub-)mm wavelengths, such 
a behavior is very unusual in the range of » 100 /xm among the 
few disks with characterized slope in this wavelength range (e.g., 
HiUentu'and et al. 2(X)8; Roccatagliata et al. 2009). Few the dkrris 
disks considered there, typical values of A at wavelengths around 
l(X)/nn are in the range of 0 to 2. Note that our Herschel obser- 
vations are the first to characterize the spectral slope of a large 
number of debris disks at wavelengths around lOO /rm. Earlier 
observations in this wavelength range suffer from low sensitiv- 
ity. Thus they were limited to few teight disks or disks with very 
flat spectral slope in the far-ir. The completion of our survey and 
an analysis of the whole data set obtained (Eiroa et al., in prep.) 
will allow to place the three disks presented in a broader context. 

We define an unusually steep SED to be an SED that exhibits 
a Arj,v 2 > 2 in at least one combination of the PACS bands of 
70pm, 100pm, and 160pm, and for which a A,,,,j = 2 would 
not be possible within the Icr uncertainties. These criteria are 
met by the three debris disks presented in this work (Table 2). 


4. SED modeling 

In this section we perform detailed analytical model fitting to the 
observed SED data to explore quantitatively the conclusions on 


the disk properties derived in the above qualitative discussion of 
the steep slope of the observed SEDs. 

4.1. Stellar properVes 

Table 3 lists some of the main stellar parameters and other rel- 
evant observational properties of HIP 103389, HIP 107350, and 
HIP 114948 (Eiroa et al., in prep.). These parameters were ob- 
tained using the DUNES Virtui Observatory tool^. The spec- 
tral types are taken from the Hipparcos catalog (Perryman et al. 
1997; van Leeuwen 2007), and the range of spectral types from 
Skiff (2010). Bolometric luminosities and stellar radii have been 
estimated firom the absolute magnitudes and the bolometric 
corrections using the measurements by Flower (1996); similar 
values are obtained using the bolometric correction procedure 
by Masana et al. (2006). Effective temperatures, gravities and 
metallicities are mean values of spectroscopic and photometric 
estimates from difierent works in the literature (HIP 103389 and 
HIP 114948: Allende Prieto & Lambert 1999; Gray et al. 2006; 
Holmberg et al. 2009; HIP 107350: Valenti & Fischer 2005; 
Fuhrmarai 2008). Radial velocities are taken from Maldonado 
et al. (2010) and Kharchenko et al. (2007). Rotational veloc- 
ities of the stars are taken from Holmberg et al. (2007) and 
Maitfnez-Am^iz et al. (2010). Periods are taken from Reiners 
(2006) and Messina et al. (2003). The Li 1 6707.8 A line equiva- 
lent widths are taken from Maldonado et al. (2010) in the case of 
HIP 107350 and our own estimates using reduced HARPS^ data 
for the other two stars; for these latest objects the EW(Li i) have 
been corrected from the Fel 6707 .4 A line contamination as in 
Maldonado et al. (2010). All three stars show logRJjj^ activity 
index levels (Gray et al. 2006; Martinez- Amdiz et al. 2010) in- 

^ ht^://sdc.cab.in(a-csic.eVdunes/searchfonn.jsp 

^ Based on observatioas made wittt the European Southern 
Observatory telescopes obtained from the ESO/ST-ECF Science 
Archive Facility. 
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dicadng that they are active stars (e.g., Maldonado et al. 2010). 
The observed ROSAT X-Ray luminosities, Z,x/Lboi. derived in 
the context of DUNES (Eiroa et al., in prep.) are also given. 

The masses of the stars in Thble 3 are estimated from their 
radii and gravities. We estimate the age of the stars (Ihble 3) 
from the logR^ index, and periods using the relationships as 
given by Mamajek & Hillenbtand (2008) and from the X-Ray lu- 
minosity following Garces et al. (2011). For all three stars these 
age [racers yield consistent values of less than ~ 500 Myr. We 
note that EW(Li i) values arc close to those of the Hyades en- 
velope with an age of 600 Myr. Rocha-Pinto et al. (2004) give 
an age of 370 Myr for HIP 107350 based on the log activity 
index, while isochrone-based ages are 2.3 Gyr for HIP 103389, 
with a range of 0.6 Gyr to 3.8 Gyr (Holmberg et al. 2009), 
3.1 Gyr for HIP 107350, with a range of 1.1 Gyr to 4.7 Gyr 
(Valenti & Fischer 2005), and 3.7 Gyr for HIP 1 14948, with a 
range of 2.1 Gyr to 5.2 Gyr (Holmberg et al. 2009). Given that 
the stars are located right on the main-sequence, stellar ages are 
difficult to estimate on the basis of isochrones, and that they are 
very sensitive to Tm and metallicity (Holmberg et al. 2009) we 
would favor in these cases the young values estimated from the 
activity tracers and rotational periods. This choice is also com- 
patible with the kinematics of the stars and their ascription to 
the young moving groups, less than 1 Gyr, of HIP 107350 to the 
local assodaticm (~ 150 Myr) or the Hercules-Lyra (~ 200 Myr) 
associaticm (Maldonado et al. 2010; L6pez-Santiago et al. 2(X)6), 
and HIP 1 14948 and HIP 103389 to the Castor Moving Group 
(~ 200 Myr) and the young disk population (< 1 Gyr), respec- 
tively. 

HIP 107350 hosts a T dwarf companion (HNPeg B, 
M = 22Mj) with a projected physical separation of 794 AU 
(w 43") from the host star (Luhman et al. 2007). This is suffi- 
ciently far away, so that the (xesence of this companion does 
not affect our present analysis. A very faint source (» 2mJy) is 
detected in the PACS 100 ftm image close to the predicted po- 
sition of HNPeg B. However, we do not expect HNPeg B to 
be detectable by our observations (» 0.6 m Jy at 8 /rm). Thus, we 
attribute this detection to a coincidental alignment with a back- 
ground source. 


4.2. Modeling the excess 

In this section, detailed analytical model fitting to the observed 
SED data is performed to e^cplore quantitatively the conclusions 
on the disk properties derived in the above qualitative discussion. 
This is done using the results from Reduction 1, since these re- 
sults give the steep SEDs. The fitting results are expected to dif- 
fer less significandy from standard solutions using the photomet- 
ric data from Reduction 2. For modeling of the three sources, we 
use SArd) which is part of the DUNES modeling toolbox (Lohne 
et al. 2011; Augereau et al., in prep.). A description of this tool 
can be found in the appendix. Pbr the model fitting of the in- 
frared excess SED, the above derived stellar properties are used. 
It is important to note that uncertainties of the stellar parameters 
as well as on the distance of the objects are not considered in 
our modeling of the excess. Thus, the derived uncertainties of 
the disk parameters are only the formal uncertainties from the 
fitting procedure. The uncertainties of the stellar parameters and 
distance would mostly affect the effective temperature and lu- 
minosity of the star and, thus, the heating of the dust - slightly 
increasing or decreasing its temperature. In particular, they will 
not affect significantly the slope of the derived grain size distri- 
buticn. 


For the comparison between modeled and observed measure- 
ments we use a reduced chi-squared 



Njof — N Alffec, 


(4) 


where N is the number of SED data points used, Fi, Ft, and cr,' 
are the measured flux, the modeled flux, and the uncertainty in a 
data point, is the number of free parameters, and Ndof is the 
number of degrees of freedom in our fitting. For this reduced^ 
a value of 1.0 is generally desirable in the fitting if Ndof is large, 
while in our case Njot ranges from 1 to 3 where the reduced 
is desirable to be as small as possible. 

Only selected data points are used. This results in a total of 
8 data points included in the fitting for each source. The SEDs 
of the sources can be seen in Fig. 3. The references for all pho- 
tometric data points available for the three sources can be found 
in Eiroa et al. (in prep.). The selection is done by the following 
criteria: 


- All flux measurements at wavelengths > lO/rm ate consid- 
ered, but no upper limits, since they do not give significant 
additional information in the present case. Measurements at 
these wavelengths are also included, if they are photospheric, 
because they exclude significant emission from the disk at 
these wavelengths. 

- Measurements at wavelengths < 10 ^m are not included, 
since they are at much shorter wavelengths than the shortest 
wavelen^ at which excess is detected (24 /rm . . . 32 /rm). 

- To account for the 5pfeer/IRS spectrum in a consistent way, 
an additional photometric point at A = 32/rm is extracted 
from these data following Hillenbrand et (2008) and is 
fitted along with the other photometric values. Therefore, 
10 consecutive data points centered at A = 32/rm are aver- 
aged. The uncertainty is computed by adding in quadrature 
the standard deviation of these data points and a 5% cali- 
bration uncertainty (SpitzerPRS instrument handbook v4.0; 
Tbplitz et al. 201 1). Tffis will prove to be sufficient to get the 
whole SpitzerfiRS spectrum properly reproduced by all the 
best-fit models. 


The capabilities of SAtiD are used to explore a very broad 
range of parameters. The following model is employed: 

- Radial surface density distribution X (r) oc r~° with iimer and 
outer cut-off radii Tu, and 

- Differential grain size distribution dn(a) oc a^^da with lower 
and upper cut-off size a„,jn and Omax, 

- Two possible grain compositions — pure astronomical sili- 
cate (Draine 2(K)3) and a 1:1 mixture of astronomical sili- 
cate and ice (Lohne et al. 201 1; Augereau et al., in prep.) - 
to explore the possibility that water ice might be a significant 
constituent of debris disk dust and might be responsible for 
the peculiar shape of the SEDs. 

Since the number of free parameters in the fitting is close 
to the number of data points used, the fitting is expected to be 
very degenerate. To explore the parameter space in an efficient 
way and to find parameters that can be fixed, because they have a 
unique best-fit result or they have no significant effect on the fit- 
ting at all, we use different approaches. These approaches are de- 
scribed below. The parameter space explored for each approach 
is listed in Table 4. The fitting results are compiled in Table 5. 
Simulated SEDs from our best-fit models are shown in Fig. 3. 
The parameter space explored and the parameters fixed in each 
approach are motivated by the results of the previous approaches 
as described in the following. 
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Table 4. Explored parameter space for the fitting based on simulated thermal annealing (for details see Sect. 4.2) 


Approach 1 Approach 2 Approach 3 Approach 4 


ParametH' 

Range 

# values 

Range 

rnmim 

Range 

# values 

Range 

# values 

distribution 

[AU] 


817 

3-100 

817 


817 

rnrao 

817 

temp 

/fc«[AU] 

5-300 

671 

5-300 

671 

5-300 

671 

5-300 

671 

temp 

a 

-2.0 -3.5 

56 

fixed 

1 

fixed 

1 

fixed 

1 

liitear 

fliriD [A®] 

0.2-20.0 

448 

0.2-20.0 

448 

0.2-20.0 

448 

0.2-20.0 

448 

logarithmic 

[pm] 

fixed 

1 

fixed 

1 

fixed 

1 

1.0-1000.0 

674 

logarithmic 

r 

2.0 -5.0 

31 

2.0-10.0 

81 

fixed 

1 

fixed 

1 

linear 

Af<hul [Af«] 

free 


free 


free 


free 

• . . 

continuous 

Conposition 0%, 50% ice 

2 

0%, 50% ice 

2 

0%, 50% ice 

2 

0%, 50% ice 

2 



4.2.1. Approach 1 : Exploring the parameter space 

In a first approach, a range of &ee parameters is considered 
that includes those used to explain most of the known debris 
disks. An upper grain size of 1 mm is chosen, large enough to 
consider any effect of large grains on the fluxes at all wave- 
lengths observed. Seven free parameters are used (r^, rout, o, 
cL^a, y. A/dust, and dust chemic^ composition). This is the largest 
number possible considering a total of 8 SED measurements. 
Although strong degeneracies are expected, this approach is used 
to explore the parameter space without strong initial ccmstraints 
that might bias the results. Parameters that have a unique solu- 
tion or that have no significant effect on the fit can be fixed in 
subsequent approaches. 

A very narrow ring structure is found to be the best-fit for all 
three dislra, although the parameters describing the spatial dust 
distribution are not constrained very well. Furthermore, a very 
steep grain size distribution with a value for the exponent y very 
close to the edge of the explored parameter space (5.0) is found 
as best-fit for all three SEDs (Table 5). In the following, it is re- 
ferred to this as an underabundance of large grains. Since this 
is an atypical, but not unexpected result (Sect. 3), more effort is 
put on the evaluation of the significance of this result in the sub- 
sequent approaches. As the radial extent of the disks has been 
found to be narrow, the exponent of the radial surface density 
distribution is fixed in the further fits to a = 0.0 (constant sur- 
face density), decreasing the explored parameter space by one 
dimension. In a narrow ring, this will not have any effect on the 
SED. 

4.2.2. Approach 2; A larger range of possible values for y 

Now, the aim is to obtain a best-fit value for y that is included 
in the considered parameter space. A range of possible values 
for y of 2.0 to 10.0 is explored. All other parameters (beside a 
which is now fixed as described above) have the same ranges as 
in Approach I . This results in a total of 6 free parameters. 

For HIP 107350, the range of explored values of y seems to 
be still too small. This is ignored, since such a very large value 
is not significant (confidence levels: 5.9 ... 10.0, Table 5). For 
the other two objects, the value of y is well within the explored 
range. The confidence levels of the inner and outer radius sug- 
gest that the ring-like shape of the disk is not very significant, 
even for a constant surface density. 

4.2.3. Approach 3: Fixing y to 3.5 

Fixing the value of y to 3.5 and having a fixed upper grain size 
of Oniax = 1-0 mm, one can force a fit where the parameters of 


the grain size distribution are consistent with an equilibrium col- 
li sional cascade (Dohnanyi 1969). The explored ranges of all 
other parameters are unchanged compared to Approach 2. This 
results in a total of 5 free parameters. 

The resulting is much worse (by a fiictor of 2.5 to 12.6, 
Thble 5) compared to the results from Approach 2. The changes 
of the radial position of the dust ring are in line with the expecta- 
tions from the much lower abundance of small (warm) particles. 

4.2.4. Approach 4: Fixed value of y = 3.5, but free upper 
grain size 

Another possibility to produce an underabundance of large 
grains in Ae model is to let the upper grain size be a free param- 
eter. Since a very steep size distribution and a small upper grain 
size have comparable results (removing large grains from the 
model), one can fix the value of y to 3.5 to be consistent with an 
equilibrium colhsional cascade. The resulting upp>er grain size is 
expected to be sufliciently small. This results in a total of 6 fi-ee 
parameters. 

This approach gives in general the best^ (Thble 5). Upper 
grain sizes of few tens of micron are found. The low radial ex- 
tend of the debris ring is again not very significEuit, but still the 
best-fit. 


4.2.5. Grid search and inclusion of spatial information for 
HIP 114948 

From our discussion in Sect 3 we are not able to constrain our 
parameter space by any reasonable physical assumptions. Since 
SAnD is designed to explore a large, high-dimensional parameter 
space using a statistical approach, it was particularly useful to 
find the global best-fit to Ae data in the above approaches. The 
over-all result of our fitting (a lack of large grains) is clearly vis- 
ible for HIP 103389 and HIP 107350. For HIP 1 14948 the result 
is not that stringent. The value of y = 4.7;|Jj from Approach 2 of 
the SAnD fitting is significantly smaller thw the values derived 
from the other two disks (although still significantly larger than 
the value of 3.5 expected fiom a standard equilibrium colhsional 
cascade, Dohnanyi 1969). Thus, we perform a grid search of the 
parameter space using GRaTer (Augeteau et al. 1999; Lebreton 
et al. 201 1) to confirm our results from the statistical approach 
of SAitD. Furthermore, the constraint that the source is spatially 
unresolved is included in this fitting. This is not expected to 
significantly change the best-fit results (since these results are 
consistent with an unresolved disk) but to further constrain the 
range of values possible within the uncertainties of the derived 
values (since very extended configurations that would clearly be 
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Fig. 3. Observed and modeled SEDs of the three sources. The SpitzerjlSS spectrum is plotted in red. The data points considered for 
the fitting are shown in dark blue. The PACS data considered are horn Reduction 1, while the results from Reduction 2 are plotted 
in light blue for comparison. Top-. Final models from the SAnD approaches. The modeled SEDs (solid line) for the three disks are 
computed from the results of Approach 4 using pure astronomical silicate (Sect. S). The dashed line in each panel represents the 
model stellar photosphere. Bottom: Comparison between the best-fit results from the different approaches. The model^ SEDs are 
computed from the best-fit result of each approach using pure astronomical silicate. The results itom Approach 1 are not included, 
since the best-fit results are found not to be in the explored range of parameters in this approach. 


resolved are also included in the confidmice levels from pure 
SED fitting). The results from fitting with and without the spa- 
tial constraints are compared. 

Since the general result from the SAnD fitting (steep grain 
size distribution, large lower grain size) is independent from 
the approach used, we decide to use the standard approach of 
GRaTer. This is a power-law approach (dn (a) « a~^da) with up- 
per and lower limit for the differential grain size distribution and 
a two-power-law approach for the radial surface density distri- 
bution: 

S(R)« yj2{(R/RoT^- + (R/Ray^y\ (5) 

where a and 0 ^ are slopes of the radial surface density distribu- 
tion outside and inside a peak position. We fix ain = 10 to get a 
sharp inner edge close to Rq comparable to the single power-law 
approach used for the SAnD fitting. 

The spatial constraints are included in our fitting using a ra- 
dial brightness profile derived by radial averaging of the sources 
image in our 100/rm data After a Gaussian is fitted to the image 
to derive the source center, the average flux and distance from 
the center is derived for all pixels in one pixel wide radial bins. 
The uncertainty on the flux in each bin is derived using the stan- 
dard deviation of the pixel values in each radial bin. Images are 


simulated from the models, convolved with the PACS PSF mea- 
sured on the bright standard star a Boo and profiles are extracted. 
These profiles are compared to the one derived from the data. 
The same number of SED data points as before is considered. In 
addition, a total of 10 profile points are used. An over-all reduced 
including SED and profile data is computed following Eq. A. 1 
(using only one profile). The explored parameter space and best- 
fit values for the fit with and without the profile included are 
given in Table 6 . These results should be compared to the results 
from Approach 2 of the SAnD fitting, since these two approaches 
are the most similar. They are found to be consistent. Maps of the 
distribution of the reduced in the parameter space for the fits 

including and not including the proWe have been produced as 
cuts through the parameter space in three planes. The remaining 
parameters are fixed to their best-fit values, respectively. These 
maps are displayed in Fig. 4. We find that Ae best-fit values do 
not change significantly when including Ae spatial information. 
This confirms Aat the results from Ae pure SED fitting are fully 
consistent wiA a spatially unresolved disk. Pig. 4 illustrates a 
number of results from Ae mapping of Ae parameter space: 

— There is a correlation between Ae lower dust grain size 
Omin and Ae exponent 7 of Ae size Astribution. However, 
from our SED fitting one can rule out Aat any combination 
of significantly smaller grains and smaller exponent of Ae 
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Table 6. Fitting result for HIP 1 14948 using the grid search method with GRaTsr 


Parameter 

Range 

Spacing 

# values 

Fit w/o spatial constraints 

Fit w/ spatial constraints 


1.0-250.0 

logarithmic 

60 

13.7 

13.7 

a 

0.0-10.0 

linear 

11 

8.0 

9.0 

amii [^ni] 

0.05 - 43.5 

logarithmic 

55 

14.1 

14.1 

r 

2.7 -5.9 

linear 

17 

4.7 

4.7 

Vte/(V„ + V,„) 

0.0 -0.9 

linear 

10 

0.6 

0.6 


1.0x10-^- 4.0x10-3 

logarithmic 

100 

2.36 X 10-’ 

2.20 X 10-* 

reduced 

- 

- 

- 

0.19 

0.12 



Fig. 4.^ maps for HIP 114948. Tbe maps are created as two di- 
mensional cuts through the^ distribution in the parameter space 
searched with GRaTer without (fop) and with (bottom) the spa- 
tial infcnmation included. The parameters not plotted are fixed 
to their best-fit value, respectively (Thble 6). The correlations 
between the three most relevant parameters (y, Omm. and Rq) ate 
shown. The white cross denotes the position of the best-fit {low- 
est The dashed vertical line denotes the value of y = 3.5 
expected from a standard equilibrium collisional cascade. 


size distribution results in reasonable fits to the data (which 
w ould be closer to the expectations from a collisional cas- 
cade and grains close to the blow-out size of the system of 
a < 1.0/rra). 

- There is nearly no correlation between the radial position of 
the dust belt and the exponent of the grain size distribution 
in both the case without and with the spatial information in- 
cluded. 

- The degeneracy between the lower dust grain size and the 
radial position of the dust belt known from SED fitting is 
dearly visible. It is also obvious that the inclusion of the 
spatial information in part breaks this degeneracy since it ex- 
cludes a large (> 70 AU) radial distance of the belt from the 
central star. 

- In particular, the results of 0^11 and y seem to be less con- 
strained when including the spatial information. This is in 


part because the spatial information gives no constraints on 
the parameters of the grain size distribution (beside the partly 
broken degeneracy of the lower dust grain size and the radial 
position of the dust which results in an indirect constraint on 
Omin). Thus, including the radial profile we include a num- 
ber of measurements that are not sensitive to changes in y, 
which results in a less significant increase of the reduced^ 
with changes of y in this case. In addition, the data points in 
the profile can not be considered to be independent from each 
other. Although this is mitigated in part by the weighting in 
Eqs. A.l to A.3, the combination of;^^ drived from SED 
and profile is not fully consistent, in particular in the case of 
unresolved sources. To measure the quantity and uncertainty 
of a model parameter through model fitting using reduced 
;t^, one should in general only include measurements in the 
fitting process that are sensitive to the parameter (as far as 
one can tell in advance) and independent from each other. 
Thus, the uncertainties on y derived firom pure SED fitting 
should be taken more reliable here. 

5. Results 

The general results and conclusions that can be drawn from the 
fitting approaches are summarized in the following. 

- Values of y = 4.7 to 10.0 or a small upper grain size are 
found rather than the expected value of y 3.S and a large 
Onux- Models wittt y = 3.5 do not reproduce the obsoved 
SEDs in a reasonable way. This result is not strongly affected 
by any modeling degeneracies found (that are included in the 
estimate of the uncertainties). 

- Evidence is found for a large lower grain size of the 
dust compared to the expected blow-out size of the sys- 
tems (0.5 /rm to 0.9/um). However, the lower grain size 
depend very much on the approach used. Earlier studies 
of debris disk SEDs had similar results (e.g., HD 107146; 
Roccatagliata et al. 2009) that have been at least mitigated 
and partially attributed to modeling degeneracies in pure 
SED fitting after including resolved data in the fitting pro- 
cess (Ertel et al. 2011). Thus, this result has to be treated 
with caution. 

- The dust in all three disks must be located at a distance from 
the star of no more than a few tens of AU as can be con- 
cluded from the modeling results and the fact that the disks 
are spatially unresolved in the PACS images. Stronger con- 
straints cm this value are not possible due to degeneracies in 
the mcxleting. 

- The disks appear to be narrow rings. However, the uncertain- 
ties - in particular of the outer disk radius - are very large 
and also very broad disks can result in fits on the SED that 
fitU within the derived 3 cr confidence levels. This can only 
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partly be excluded by the additional constraint that the disks 
are spatially unresolved. 

It is important to note that these results depend in part on the 
model used and on the parameter space explored. For example, 
the composition and shape of the dust grains have not been ex- 
plored in detail in our modeling, but might allow one to repro- 
duce the SEDs with models that are more in line with a standard 
equilibrium collisional cascade. 

Furthermore, the results from Reduction 2 have not been 
considered in the fitting. Considering these photometric mea- 
surements and uncertainties instead of those from Reduction 1 
would lower the significance of the underabundance of large 
grains. As a result, values are expected to be found for the model 
parameters that are more in line with those expected from an 
equilibrium collisional cascade. A^thout a deeper understanding 
of the data obtained with Herschel, it is not possible to make a 
final statement about the significance of our modeling results. 

Each approach presented results in a different set of best-fit 
parameters and a different;^, while the results can be interpreted 
in a consistent way. To find from the fitting results a final model 
of each disk, the following selection is ^phed: 

- The results fix)m Approach 1 have to be ruled out, because 
the best-fit parameters from this approach are close to the 
edge of the explored parameter space. 

- In general, the SBDs cannot be fitted by a grain size dis- 
tribution that includes a significant amount of large grains 
(Approach 3). 

- The results from Approach 2 (free y, an.„ = 1.0 mm) and 
Approach 4 (free Omax, y = 3.5) are consistent, while the re- 
sults from Approach 4 and pure silicate represent in general 
the best fit to the data. Thus, these results are considered as 
the best-fit models from the fitting. 

From these final models, one finds dominating dust temperatures 
(the temperature of the smallest grains that dominate fte emis- 
sion) of63K, 50 K, and 73 K for HIP 103389, HIP 107350, and 
HIP 114948. 

5. 1. The origin of the dust 

The existence of micron-sized dust grains around main sequence 
stars is usually explained by collisions of larger bodies such as 
planetesimals (e.g., Krivov 2010). In such a scenario one would 
expect a significant amount of larger grains (several tens of mi- 
cron to millimeter sized particles) to be present that are pro- 
duced by the collisional cascade also producing the smallest 
grains (Dohnanyi 1969). These grains significantly contribute to 
the SED in the observed wavelraigth regime. Assuming that the 
steep slopes of the SEDs are real, we find from our modeling 
strong evidence for an underabundance of these larger grains. 
To our knowledge, these are the first debris disks discovm'ed, 
that potentially exhibit such a peculiar shape of the SED in the 
wavelength regime of 70fim to 160^m. All this would sug- 
gest that these objects are exceptions from the common under- 
standing of dust creation in debris disks. On the other hand, six 
more disks vrith similar shape of the SED have recently been 
identified in our DUNES survey, which might suggest that this 
phenomenon is common among low-mass (i.e., low luminosity, 
Thbles 5 and 7) debris disks^. The fact that they arc very faint 
(Thble 7) might imply that we are faced with a new class of de- 
bris disks that were not revealed earlier due to limited sensitivity 

’ For an oveniew of the DUNES results we refer to Eiroa et al., in 
prep, 


in the relevant wavelength range. In the following, we discuss 
scenarios that might be capable to explain the presence of a disk 
as it has been modeled for the three systems. 

Scenario 1: Significant deviation from the conditions re- 
quired for a standard equilibrium collisional cascade 
For the standard equilibrium collisional cascade, a number of as- 
sumptions are made that are not necessarily valid in debris disks 
due to the effects of radiation pressiue and Poynting-Robertson 
(PR) drag (Sect. 1). For the massive debris disks known so far, 
PR drag is expect^ to be negligible (Wyatt 2005). Following 
Backman & Paresce (1993), one can estimate the ratio between 
colUsional time scale and PR time scale for the three disks 
to 10“^ . . . 10“^. However, the disk models used (best-fit from 
Approach 4, pure silicate) represent very narrow rings, which 
is not a significant fitting result. Assuming a ring with a width 
of 10 AU starting at the inner radius found from the fitting, one 
finds tcou/tpR = 3 X 10“^ . . .5 X 10"^. It is not possible to put 
strong constraints on the dust dynamics due to the crude esti- 
mate of the time scales. Thus, transport mechanisms may play a 
role in the dust dynamics of the small grains in these disks. The 
effects of these processes have been modeled, e.g., by Krivov 
et al. (2006); Thebault & Augereau (2007); Wyatt et al. (201 1). 
In particular, radiation pressure leads to a wavy structure of the 
grain size distribution increasing the slope of the grain size dis- 
tribution in certain ranges . PR drag results in a depletion of small 
grains, increasing the dominating grain size, but mitigating the 
wavy structure of the size distribution. 

Scenario 2: Different chemical composition or physical shape 
of the dust grains than the assumed one 

One might imagine a dust composition emitting significantly 
more efficiently than astronomical silicate or ice in the wave- 
length range of 70fim to lOO^rm compared to longer wave- 
lengths. Therefore, one would need grains with a break in 
Gabs W (Sect. 3) at shorter wavelengths (in particular for large 
grains). Also, if the large grains were significantly colder than 
expected, their emission would be reduced. This would both al- 
low a significantly larger amount of large grains to be present in 
the systems than modeled. Furthermore, Mie calculations might 
result in an incorrect approximation of j2»bs particularly for large 
grains (tens of microns or larger), e.g., due to incorrect assump- 
tions on the physical shape of the grains. 

Voshchiimikov et al. (2006) simulated the effect of poros- 
ity on the absorption efficiency of dust grains. They found that 
the temperature of dust grains at a given distance from a star 
decreases significantly with increasing porosity. If one now has 
small, compact grains and larger, porous grains, the difference 
in tempm'atures of the smaller and larger grains can be increased 
significantly. This would be the case, if the large grains produc- 
ing the smdler ones through collisions were porous, composed 
of smaller, compact units in the order of the lower grain size 
derived from the modeling of the three systems (i.e., 5/rm to 
10^). Such a scenario would result in an overabundance of 
grains of this size compared to the larger, weaker grains. Such 
an overabundance would shift the lower cut-off of our adopted 
power-law distribution toward this grain size, although smaller 
grains are present in the system. 

Scenario 3: A shepherding planet 

Planets can trap dust particles into mean-motion resonances 
(e.g., Wyatt 2006). This results in a barrier against particles mov- 
ing inward due to PR drag. This barrier is less efficient for very 
small grains for which PR drag is very strong (Reidemeister et al. 
201 1). On the other hand, very large grains and planetesimals are 
not significantly affected by PR drag. 
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Ibble 7. Simulated observational properties of the disks derived from our final models 


Source Aim/L> 





du»l / ^ 

fac&oD 

edge-on 

iBce-on 

edge-on 

face-cm 

edge-on 

face-on 

edge-on 

HIP 103389 1.5xl(r= 

1.1 X 10“* 

1.9x10-* 

1.1 X 10““ 

1.7 X iir* 

1.3 xur* 

1.8 X 10“* 

1.8 X 10-* 

1.5 X 10-* 

HIP 107350 0.6 X 10“’ 

4.9 X 10"’ 

7.4 X 10“® 

4.6 X 10“’ 

8.5 X 10“« 

5,5 X 10“’ 

1.3x10“’ 

6.9 X 10“’ 

3.9 X 10“’ 

HIP 114948 2.5x10“’ 

2.1 X 10“’ 

3.4 X 10"* 

2.0 X 10“« 

3.0 X 10-* 

2.6 X 10-*> 

3.3 X 10-* 

3.7 X 10-® 

2.6 X 10“* 


The dust seen in these disks might be produced in a faint, 
transport dominated (Krivov 2010) debris disk further away 
from the star, too cold and too faint to be detected in the available 
data. The small particles (few tens of micron and smaller) would 
then be dragged inward by PR drag. A possible planet present 
further inward in the system should trap the particles into reso- 
nance. This would result in an accumulation of particles with a 
very distinct range of sizes, which would explain both the lack 
of large grains (at least in an abundance and at a position where 
they wotild be detectable) and the lack of small grains. Such a 
“dust trail” has been observed, e.g., to be associated with our 
Earth and has been modeled to predominantly consist of grains 
of w 12/rm and larger (Dermott al. 1994). It will have to be 
evaluated, whether an outer planetesimal belt capable to produce 
enough dust through collisions to replenish the dust in the disk 
can be faint enough not to be detected by the Herschel obser- 
vatiens, and under which conditions the planet can trap enough 
dust into resonance to produce the disk observed. This has to be 
done through detailed dynamical and collisional modeling (e.g., 
Reidemeister et al. 2011), which is not within the scope of the 
present work. 

5.2. Observational perspectives 

Further observations with present and near future instruments 
can help to reasonably increase our understanding of these debris 
disks. In the following, we discuss a number of observations that 
might result in very valuable data. 

- Since the most extended models of the three disks predict 
radial extents of only w 2'.'2, these objects arc at the edge 
of what is resolvable with coronagmpky. Furthermore, the 
disks are very faint in scattered light (l^le 7). It is impor- 
tant to note that in the case of edge-on orientation most of the 
Gux comes from forward scattering and will then be concen- 
trated close to the star (projected physical separation) with 
no contribution to the signal in coronagraphic observations. 
A contrast ratio of < 10“® is not accessible to present instru- 
ments. 

- Successful opticai/near-infrared imaging of planetary com- 
panions would give strong evidence that Scenario 5 (shep- 
herding planet) is responsible for the peculiarities of the disk. 
Determining the position of this planet would also help to 
further constrain the position of tte dust in this scenario. The 
youth of the stars means good chances to directly image gi- 
ant planets at separations of > 1" from the star with present 
methods (e.g., Marois et al. 2010). 

- The predicted extents of the disks are in a range easily resolv- 
able with ALMA. However, the disks are already very Mnt 
at PACS wavelengths and the expected surface brightness 
is decreasing particularly steep toward wavelengths acces- 
sible with ALMA. It is not clear, whether ALMA observa- 
tions ate sensitive enough to detect these disks. On the other 


hand, at (sub-)mm wavelengths even upper limits in the sen- 
sitivity range reachable are expected to provide useful fur- 
ther constraints on the shape of the SEDs. From the results 
of our fitting Approach 4, predicted fluxes of star and disk 
at 350 /rm (the shortest wave band accessible with ALMA) 
are 0.55 mJy and 0.11 mJy for HIP 103389, 0.53 mJy and 
0.05 mJy for HIP 107350, and 0.60 mJy and 0.20 mJy for 
HIP 1 14948, respectively. 

- Further photometry and spectroscopy in thermal reemission 
would help to constrain the radial distribution of the dust 
and the size distribution of the small grains (wavelengths 
between 40 pm and 70 pm) and the steep slope and the 
shape of the SEDs at longer wavelengths (A > 70 pm). The 
Stratospheric Observatory for Infrared Astronomy (SOETA) 
is most promising to provide observational capabilities in the 
relevant wavelength regime in the near future. 

6. Conclusions 

The first data at wavelengths > 70 pm for the three debris disks 
discussed in this paper have been presented and modeled. All 
three sources potentially exhibit an unusually steep decrease of 
the SED in the wavelength range between 70pm and 160^m. In 
a general discussion, it has been shown that this peculiar shape 
of the SED is an indicator for a deviation from the case of a stan- 
dard equilibrium collisicmal cascade assuming standard grain 
composititxi (astronomical silicate or ice) and shape (compact, 
spherical grains). However, it has also been shown that the un- 
derstanding of faint source photometry obtained with Herschel 
is still incomplete and that the results presented here depend 
very much on the actual version of the data reduction pipeline. 
Provided that the steep decrease of the SEDs is real, modeling 
implies that the therm^ emission from these disks is dominated 
by a very distinct grain size regime of several micron to few 
tens of micron. The disks have been modeled as narrow rings 
with a significant underabundance of large grains. A number of 
possible explanations for such an unusually steep shape of the 
SEDs have been discussed. Six more candidates for this new 
class of debris disks have been identified so far from the on- 
going DUNES survey. This is the first published discovery of 
debris disks that potentially exhibit such a peculiar shape of the 
SED in this wavelength regime. 
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Appendix A: Description of SAnD 


When performing analytical model fitting of debris disks, one 
is faced with a number of challenges. On the one hand, re- 
cent results (e.g„ Backman et al. 2009; Ertel et al. 201 1) reveal 
more and more the complexity of the known debris disks and 
the inability to reproduce the available data with simple mod- 
els. Furthermore, standard assumptions (e.g., lower grain size 
consistent with the blow-out size of the system, grain size dis- 
tribution exponent of -3.5) have been proven to be inadequate 
(e.g., Ertel et al. 2011; this work), and unexpected results like 
an outward increasing surface density distribution (Ldhne et al. 
2011) have been found. The increasing number and quality of 
the data available (in particular, spatially resolved images tom 
Herschel and ALMA) allow one to break modeling degenera- 
cies through simultaneous multi-wavelength modeling and to in- 
crease the complexity of die models. On the other hand, finding 
a best-fit model in a complex parameter space results in signif- 
icant challenges on the fitting method and the simulation of the 
observed data in order to reduce the computational effort neces- 
sary. 

Tb rise to these challenges we, developed the tool SAnD. It is 
able to fit SED data and radial profiles in thermal reemission si- 
multaneously using a simulated annealing approach (Press et al. 
1992) on a grid of possible values in the parameter space. It uses 
an analytical approach for the radial density distribution of a ro- 
tationaUy symmetrical disk as well as for the grain size distri- 
bution. The SED is computed in an analogous way to Wolf & 
Hillenbrand (2005). Resolved images are computed with suffi- 
cient resolution and convolved with the telescope PSF. Radial 
profiles are then extracted along the major and/or minor axis. 
Each image (i.e., each pair of major and minor axis profiles or 
single majm: axis profile, if no minor axis profile is extracted, 
e.g., in azimuthally symmetrical disk images) is scaled by a fac- 
tor Xi to minimize the;^f^ (Eq. A.3). This way, only the shape 
of a radial profile is fitt^. This is done because flux calibration 
uncertainties have to be considered only oiKe per wavelength, 
but would be included for each data point in the radial profiles 
if absolute profiles were considered. I^ithermoie, deficits of our 
models to reproduce the absolute flux at one wavelength at which 
also radial profiles are fitted (e.g., due to uncertainties in the op- 
tical properties of the dust) would be considered in each profile 
data point and the SED data point which would result in an over- 
weighting of the data at this wavelength. Simulated and observed 
multi-wavelength data are then compared using the following re- 
duced 




with 


and 



M 


(A.1) 

JVu+^uJ’ 



(A.2) 




“) ■ 


(A.3) 


where « = 0 is for the SED and « > 1 for the profiles, / = 1 for 
major axis and j = 2 for minor axis. In equations A.l to A.3, the 
notations used are; 
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- M - number of images 

- Nq- number of SED data points 

- iV(,i - number of profile points along major axis 

- - number of profile data points along minor axis 

M 2 

- Nfoi - total number of data points, Wq + 2 S Nu 

i=l^l 

“ AW - number of free parameters 

- f*. Pk - modeled and observed SED, respectively 

- Si,j,k, §i,j,k - modeled and observed surface brightness pro- 
files 

- (Tic, {Tijjc - observed uncertainties 

- Xi- surface brightness profile scaling factor 

- wi- weight. 


For the present work, only one profile is used and wo and 
vvi are set to 1.0, which means that the profile is given the same 
weight as the SED. 

In order to allow one to find the best-fit model within a given 
range of parameters of a high dimensional parameter space in 
reasonable time, one has to compute the value of ^ for a given 
set of parameters very fast. Creating the SED and three profiles 
at different wavelengAs for one particular model and comparing 
them with the obsKvaticms usually takes ~ 0.1 sec on an Intel 
Xeon E5410 CPU (2.3 GHz). The exact time depends mostly on 
the particular set of parameters and the number of wavelengths, 
at which radial profiles are provided. It hardly depends on the 
number of data points in the SED and the profiles. 

The simulated armealing approach enables us to handle a 
very large and high dimensional parameter space without sam- 
pling the whole range of parameters. Therefore, a random walk 
on the grid of possible parameters is started at an arbitrary posi- 
tion. The probability p to go a certain step depends oa;^ of the 
actual and;^^ of the the previous step. It is 1 for^ (fil be- 
comes better) and follows a Boltzmann distribution for^ <x^ 
(fit becomes worse); 



In this equation, the quantity 7 is a parameter comparable to the 
Temperature of an “annealing” physical system and controls the 
ability to reach areas in the parameter space that give worse fits. 
It is chosen at the beginning of each run to be a large value based 
on experiences from previous runs and is lowered with each suc- 
cessful step of the random walk. The width of each step is chosen 
from a probability distribution that prefers short steps, while the 
maximum step width is 1/2 of the whole range of values for each 
parameter. The probability to go long steps is lowered with each 
successful step. \^th this approach the code is able to reach each 
position in the parameter space, but prefers regions with low 
The run stops, when a position in the parameter space with the 
following properties is reached: 

- With the actual value of 7 the code was not able to leave the 
position after computing 500 further models. 

- The actual ^ is lower than or equal to the best reached 
during the whole nm. 

- The actual is lower than or equal to a maximum^ based 
on experiences from previous runs (input parameter). 

A set of model parameters that satisfies all of these conditions is 
considered as a likely global fit. To increase confidence in this 
result, one can start several runs at different regions of the pa- 
rameter space (4 in the recent case for the final fitting). Only if 


aU runs reach the same region in the parameter space, we con- 
sider the result to be really a global best-fit. If only the first and 
not more tiian one other of these stof^ing criteria is satisfied, the 
value of 7 is increased and the code goes on trying to leave the 
actual position in the parameter space. 

Error estimates are done in the context of the simulated an- 
nealing approach by starting a new random walk at the best-fit 
position with a fixed value of 7 (usually 10% of the best-fit 
in our case). The code counts how often each value of a param- 
eter is reached. This way, the code samples the projection of the 
probability distribution of the random walk (Eq. A.4) on the axis 
of each parameter. With the knowledge of the value of 7 one can 
then compute the probability distribution 

p’ (A.5) 

and corresponding levels of confidence. 

As a feature of the simulated annealing approach the time 
it takes to find a best-fit (and uncertainties) does not vary very 
much with the number of fiee parameters or the total size of the 
parameter space to be searched, but de]^nds on the shape of the 
distribution in ttie parameter space. Weak variations of the^ 
over a broad range of values of the parameters as well as very 
degenerate fitting problems require more time, i.e. more mod- 
els to be computed. For the recent work (SED fitting only) we 
find the best-fit in a parameter space composed of ~ 10'° to 10“ 
knots typically in less than 24 hours including the estimates of 
the uncertainties. Therefore the code computes about 100.000 to 
1.000.000 models and accepts (goes) about 10.000 and 100.000 
steps for the fitting and computes ~ 1.000.000 models and ac- 
cepts ~300.000 steps for the error estimates. 
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Tbble 5. Results from the fitting based on simulated thennal annealing (fen- details see Sect. 4.2) 


Parameter 




Best-fit value [3cr 

confidence levels] 




HIP 103389 


Approach I 

Approach 2 

Approach 3 

Approach 4 


Silicate 

Mixture 

Silicate 

Mixture 

Silicate 

Mixture 

Silicate 

Mixture 

Rta[AU] 

‘Ib.l [8.5- 32.2] 

20.8 [9.1-35.9] 

18.2 [7.9 -24.5] 

20.9 [8.0 -26.6] 

11.9 [4.2- 16.0] 

13.6 [4.4- 19.1] 

42.3 [12.9 -63.3] 

223 [8.6-43.0] 

/?out [AU] 

20.7 [15.3-53.1] 

20.9 [15.4-123.3] 

20.0 [16.3-65.1] 

20.9 [17.3-77.3] 

12.0 [8.8 -37.4] 

13.6 [10.2-38.5] 

46.0 [21.3 -138.5] 

22.5 [17.9-110.9] 

a 

0.1 [-2.0 -3.5] 

2.6 [-2.0 -3.5] 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 


6.6 [4.7 - 8.6] 

9.2 [4.4 -11.4] 

9.5 [7.8 -10.4] 

12.8 [10.7-13.6] 

6.1 [2.8-9.81 

6.6 [2.6- 12.0] 

4.2 [3.1 -7.8] 

9.3 [3.7 -14.9] 

rw [/rm] 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

14.3 [12.7-18.1] 

22.7 [14.3-28.5] 

y 

5.0 [4.6 -5.0] 

5.0 [4.5 -5.0] 

7.4 [6.3 -10.0] 

9.0 [6.5 - 10.0] 

3.5 (fixed) 

3.5 (fixed) 

3.5 (fixed) 

3.5 (fixed) 


1.98e-5 

1.67e-5 

1.30e-5 

1.29e-5 

4.69e-5 

4.03e-5 

4.86C-5 

1.45e-5 

reduced^ 

4.112 

3.664 

0.816 

0.776 

10.277 

8.488 

0.628 

0.760 

HIP 107350 


Approach 1 

^)proach2 

Approach 3 

A]^)roadi4 


Silicate 

Mixture 

Silicate 

Mixture 

Silicate 

Mixture 

Silicate 

Mixture 

RtalAU) 

15.3 [3.8-28.5] 

21.5 [7.4 -93.0] 

29.1 [7.9-47.2] 

30.6 [5.6-44.1] 

9.6 [3.0- 15.6] 

IM [3.()-l6.5] 

37.1 [4.3-54.4] 

35.2 [7.2 -543] 

Ro-[AU] 

15.9 [93 -227.3] 

21.8 [10.3 - 276.6] 

313 [13.9- 113.9] 

323 [16.0-138.5] 

9.6 [5.8 -33.0] 

11.0 [9.1 -35.7] 

37.4 [19.2- 1873] 

353 [17.4- 145.3] 

or 

-0.9 [-2.0 -3.5] 

1.4 [-2.0 -3.5] 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 

Omiii [;rm] 

5.9 [4.7 -9.0] 

5.9 [43 -9.1] 

6.9 [2.7 -10.9] 

8.2 [3.6- 10.7] 

5.7 [4.8-11.0] 

5.8 [4.9 -11.7] 

7.8 [1.6- 10.4] 

9.6 [2.2 -13.5] 

Odk lAtm] 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

7.8 [6.3 - 13.3] 

9.6 [6.4 -17.4] 

y 

5.0 [4.0 -5.0] 

5.0 [4.0 -5.0] 

10.0 [6.0- 10.0] 

10.0 [5.9 -10.0] 

3.5 (fixed) 

3.5 (fixed) 

3.5 (fixed) 

3.5 (fixed) 

A^dml [A^*] 

4.59e-6 

6.03e-6 

l.lle-5 

9.29e-6 

1.05&-5 

9.09e-6 

1.54e-5 

1.08e-5 

reduced 

7.048 

5.672 

1.652 

1.568 

4.229 

3.869 

1.528 

1.488 

HIP 11 4948 


Approach 1 

Approach 2 

Approach 3 

Approach 4 


Silicate 


Silicate 

Mixture 

Silicate 

Mixture 

Silicate 

Mixture 

R. [AU] 


13.6 [6.8 - 14.5] 

12.8 [7.1 - 13.7] 

13.5 [7.8- 14.9] 



323 [9.0-40.1] 


RcutlAU] 

12.8 [12.1-39.8] 

13.8 [12.9-34.7] 

12.8 [12.1-23.2] 

13.8 [12.9-24.9] 

12.7 [10.5-25.6] 

14.1 [11.0-30.8] 

34.5 [26,9 -81,4] 

13.8 [12.9 -22.9] 

a 

2.9 [-2.0 -3.5] 

0.0 [-2.0 -3.5] 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 

0.0 (fixed) 

a™, [/mi] 

9.6 [8.9 - 10.2] 

12.9 [11.9-13.8] 

9.6 [8.9 -10.2] 

12.9 [11.9 -13.8] 

6.4 [3.6 -9.9] 

7.4(3.0-11.7] 

3.2 [2.6 -4.6] 

11.8 [10.7 -13.1] 

<W [/un] 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

1000.0 (fixed) 

24.2 [10.0 -273] 

43.8 [10.0-53.2] 

y 

4.7 [4.5 -5.0] 

4.7 [4.4 -5.0] 

4.7 [4.4 -5.2] 

4.7 [4.4 -5.1] 

33 (fixed) 

3.5 (fixed) 

3.5 (fixed) 

33 (fixed) 

[W.] 

2.19C-5 

2.16e-5 

2.19e-5 

2.16C-5 

8.82e-5 

7J6e-5 

5.93e-5 

I.6IC-5 

reduced;^ 

0.568 

0.464 

0.284 

0332 

3.291 

2J33 

0.132 

0.232 


S. Brtel et al.; A peculiar class of debris disks from Herschel/DUNES 






